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Abstract 



The Plebariski-Demiariski solution is a very general ajdally symmetric analytical 
solution of the Einstein field equations generalizing the Kerr solution. This solution 
depends on seven parameters which under certain circumstances are related to mass, 
rotation, cosmological constant, NUT parameter, electric and magnetic charges, and 
acceleration. In this paper we present a general description of matter wave interferom- 
' etry in the general Plebairski-Demiairski black hole space-time. Particular emphasis 

^ ' is placed on a gauge invariant description of the symmetries of the gauge field. We 

\l , show that it is possible to have access to all parameters separately except the accel- 

eration. For neutral particles there is only access to a combination of electric and 
' magnetic charge. 

PACS: 03.75.Dg, 04.80. Cc, 04.20.Jb 

O 

O ■ 1 Introduction 

A given space-time can be explored by dynamical systems like massive point particles and 
■ light rays and, indeed, most of the interpretation of the properties of a given space-time 

H I or a particular solution of Einstein's field equation is based on that [T] . A further system 

which may be used is quantum matter. With quantum matter not only the trajectory of the 
particle is subject to observation but the amplitude and the phase of the quantum systems 
evolving in the given space-time. In its quasiclassical limit the phase is proportional to the 
proper time along the classical trajectory and, thus, contains all the information about the 
space-time metric. This phase can be accessed through interferometry. Since the phase is 
very sensitive to external fields like the gravitational field, interferometry is a very precise 
tool to explore the properties of space-times. 

The first interference experiment sensitive to the gravitational field has been carried 
through with neutrons in 1975 by Colella, Overhauser and Werner j2j. They used a macro- 
scopic interferometer: after coherent splitting neutrons moved along different paths with a 
height difference of several cm. With neutron interferometry also the rotation of the Earth 
has been observed [3]. This is a matter wave analogue of the famous Sagnac effect. Both ef- 
fects have also been observed with much higher accuracy using atomic beam interferometry 

HIS]. 

In all these cases it is sufficient to calculate the phase shift on a non-relativistic level, 
that is, by using the Schrodinger equation. However, as the accuracy of these devices 
increases and also interferometry with ultracold atoms with much higher accuracy in space 
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is feasible it might be interesting to study these interference experiments in a broader 
context. This broader context might be (i) to include relativistic effects and effects due 
to the spin of particles [H [71 [S] , (ii) to include terms violating the Einstein Equivalence 
Principle [9] , (iii) to include a general space-time metric in the frame of a post-Newtonian 
parametrized test theory [10], or (iv) to use a particular space-time metric given by a 
certain solution of Einstein's equations. In this paper we calculate the phase shift for a 
charged scalar field in a Plebahski-Demiahski space-time. 

The Plebaiiski-Demiahski family of solutions of the Einstein field equations associ- 
ated with the gravitational fields of isolated massive objects [11] are known to completely 
exhaust the Petrov type D space-times. These axially symmetric solutions are character- 
ized by seven parameters which under certain circumstances are related to mass, angular 
momentum, cosmological constant, electric and magnetic charges, NUT parameter and 
acceleration. This includes black hole space-times like Kerr-NUT-(A)dS space-times, the 
C- metric describing accelerating sources and non-expanding solutions of Kundt's class. 
See |12| for a review and a new form of this family of solutions. The non-accelerating 
Plebahski-Demiahski solutions possess the outstanding property that they allow separable 
Hamilton-Jacobi equations and, thus, integrabihty of the geodesic equation [T3l fT4l [T5]. 
This also extends to higher dimensions [HI \17\ [18] . Higher dimensional solutions of this 
type became popular in the connection of string theories and brane world models [19| . Here 
we describe charged particle interferometry in such general space-time models in order to 
answer the question whether all parameters characterizing the Plebahski-Demiahski family 
of solutions are accessible through such a type of experiment. 

A scalar field obeying the Klein-Gordon equation is the simplest quantum object. In- 
terference with quantum fields with inner degrees of freedom like spin-^ fields explore the 
same properties but are more complicated to describe. Due to their spin-curvature cou- 
pling spin-i fields are more important for the local exploration of the space-time curvature 
[20j . However, since in ordinary situations within the Solar system the space-time curva- 
ture and, thus, the spin-curvature coupling is very small, these inner degrees of freedom 
effectively do not contribute to the observed phase shift. Therefore, for standard situations 
interferometry with scalar fields is by far sufficient. In more general geometries spin-i 
fields play a distinctive role in the search for non-metrical fields like torsion This is 
not considered here. 

In this article we describe the interference of a scalar field in a Mach-Zehnder type 
interferometer. The outline of the paper is as follows: First we solve the Klein-Gordon 
equation in a general Riemannian space-time within a quasiclassical approximation. Then 
we present the mathematical description of the Mach-Zehnder interferometer which has 
been realized by neutron as well as atomic interferometry. In relativistic terms such a 
Mach-Zehnder interferometer forms a worldtube, see Fig. [T] The combination of the prop- 
agating wave and the interferometer geometry results in the observable phase shift. This 
phase shift is then specialized to a Plebahski-Demiahski space-time which depends on seven 
parameters. We discuss how the various parameters of the metric contribute to the cal- 
culated phase shift. This represents a generalization of the description of interferometery 
of neutral particles in this class of space-times t21j. In our approach we put particular 
emphasis on a gauge invariant formalism for the description of space-time symmetries of 
the electromagnetic field. This formalism naturally leads to, e.g., gravitationally modified 
gauge invariant electrostatic and magnetostatic potentials. These U{1) gauge invariant 
potentials in turn define a restricted gauge freedom which is uniquely related to the choice 
of a scale and a zero of the energy. This describes the fact that experimentally only ratios 
of differences of energies can be measured. 

2 Dynamics of the matter field 

We take as quantum system a scalar field with electric and magnetic charge. This scalar 
field is assumed to obey the Klein-Gordon equation minimally coupled to the electromag- 
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Figure 1: A rotating and accelerating Mach-Zelinder interferometer defines a worldtube. Each constituent 
of tlie interferometer, e.g., tlie beam splitter, the mirrors, the analyzer, moves on a 4— worldline in In a 
stationary space-time each of these worldlines are one point in the quotient space yl^ which in this Figure 
is isomorphic to the x — j/-plane. 



netic field and to gravitation 

g^"" - m^iP =^ . (1) 

Here g^i, is the non-singular space-time metric with signature +2 and g^'^ its inverse 
defined by g^''gpu — S^. We also set ?i = c = 1; greek indices run from to 3. The 
covariant derivative % is defined by %T'' = D^T'' - iepA^T" - gpA^T", with D^T" = 
dfiT" + {^p} TP where {Ifp} = ^g^"" {d^gpa + dpg^cr - d^g^'p) is the ChristofFel symbol and 
Cp and gp are the electric and magnetic charges of the field T^. The Maxwell potential Ap 
gives the Maxwell field strength = dpA,j — d^Ap. In vacuum dpF'^'^ = and dpF^^" — 
where F^" := j^e'^'^P'^ Fp„ is the dual field strength and e^'''"^ the Levi-Civita symbol. 

The dual potential Ap is related to the dual field strength, see e.g. [22] ■ Note that A and 
F are purely imaginary. 

The solution of the Klein-Gordon equation we are looking for can be obtained through 
the ansatz ipix) — a(x)e*'^'-^-' where S is assumed to be real valued and a to be invariant 
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under gauge transformations. We insert this ansatz into the field equation and obtain 

-a {d^S - epA^ + igpA^) {d^S - CpA^ + igpA^,) e^^] - m^ae'^ (2) 

With = d^S and = — CpAf^ + igpAf^ this reads 

= g^''D^D,a + tg^-'{2D,ap^ + D^p,a)-gt'y,p^a~m^a. (3) 

Now we assume that the amplitude a and the kinetic momentum vary very slowly 
compared to the variation of the phase S. That means 

df,a < p^a , g'^'D^p^ < . (4) 

These are the conditions for the quasiclassical approximation. If we take h ^ 1 then these 
conditions can be mimicked by the requirement that h is small. We will come back to these 
approximation conditions. 

With ^ we obtain from ([3]) to lowest order 

= (-.g^>.PM - m^) a , (5) 

that is, 

g'^^'Pt^P,^ = -"1^ . (6) 

This is the general relativistic Hamilton-Jacobi equation for a point particle with mass m. 
Therefore p^ and can be interpreted as kinetic and canonical momentum, respectively. 
From the one-form we obtain the 4-velocity vector through = -^g^^Pv The integral 
curves of v are realized by the trajectories of the peaks of wave packets. Differentiation of 
(pl gives the Lorentz-force equation for a particle with electric and magnetic charge 

mDyV — epF{v) — igpF{v) . (7) 

If ([6|) is fulfilled then we obtain from ([3]) 

= gf^'D^D^a - i (2g'^%i?,a + g^^'D^p.a) . (8) 

Owing to the higher derivative the first term is smaller than the rest so that we obtain as 
first order approximation 

= 2g^''p^D,a + g'^'D^p.a . (9) 
We replace the momentum by the velocity and obtain 

D,a^-^ea, (10) 

where 6 = Df^v'^ which is interpreted as expansion of the congruence of worldlines given by 
the integral curves of v. This is a propagation equation for the amplitude a. This equation 
can be rewritten as 

— — a=— lna = (11) 

a as as 2 

with the solution 

a = exp y ()ds^ ag (12) 

where ag = a(a;o) with the initial point xq. 

The quasiclassical wave function then is given by 

ip{x) = exp ( - ^ / Ods] exp ( -i [ S^^dx^"] a^ (13) 
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where the integration is along the trajectories given by the solutions of ([7]). 

Now we come back to the approximation condition ([?]). The second condition im- 
plies 9 <^ m — 1/Ac, where Ac is the Compton wavelength. In other words, \cO <C 1. 
The expansion 6 characterizes the change of the spatial distance between two neighboring 
trajectories: if 5r denotes a spacelike vector between two neighboring trajectories then 
Dtjlnp ~ ^0 where = g{6r,Sr). As a consequence our condition reads XcD^lnp <C 1 
and means that the relative change of the distance between neighboring trajectories is very 
small within a time span of a Compton time. Using the Raychaudhury relation this also 
implies a condition on the strength of the space-time curvature and the electromagnetic 
field. 



3 Model of the interferometer 

3.1 Stationary congruences 

A congruence of Killing trajectories is a set of worldlines whose tangent vector is propor- 
tional to a Killing vector field. A Killing field ^ fulfills = where is the Lie 
derivative with respect to ^. This is equivalent to Df^(_„ + D^S^^ = 0. 
If M is a normalized vector field proportional to ^ we have 

C = e%, (14) 

where g{u,u) = —1. The exponent [/ is a generalized gravitational potential since its 
gradient gives the acceleration of the Killing trajectory, a — DuU — g{-,dU). It is easy 
to show that this potential is constant along the Killing trajectories, DuU = 0. The 
generalized potential U combines gravitational as well as inertial forces. 
We also need the projector defined through 

P,A = A^'-^^, (15) 

which projects a vector A onto the rest space of u. It is — and — Pu. 

The rotation u^i, of a congruence is defined by u^i, — {Pu)'^{Pu)1D[pUcr]. In terms of 
the Killing field we obtain [23] 

i^fiiy = -e^^-Df^Cn + 29[^f7u^] . (16) 

Note that lo^^S^'^ = and ^^uj — 0. The stationarity of a congruence is equivalent to three 
conditions which have a direct physical interpretation: (i) the congruence is rigid, (ii) the 
rotation vector in the rest frame does not precess, and (iii) the acceleration in the rest 
frame of the congruence is co-rotating [25] . 

The existence of Killing trajectories in a 4-dimcnsional manifold ^ allows to establish 
a 3-dimensional quotient space jV with an induced positive definite metric h^^ [24[ 125] . 
This metric on -jV and its inverse are given by 

In adapted coordinates the covariant form of the metric tensor adapts the form hij = 

9ij - goigoj/goo [26]. 

3.2 Stationarity of the electromagnetic field 

In addition to the stationarity of the interferometer we also require the stationarity of the 
electromagnetic field given by = 0. The Lie derivative acting on differential forms 

employs the identity — i^d + di^, where d is the exterior derivative operator and is 
the interior product. With the homogeneous Maxwell equations Jif^F = this implies that 
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££^A = d\ for some function A (here one uses = d££^). Since for Killing vectors 
commutes with the duality operation we also have = 0. With the Maxwell equations 
in vacuum this yields — dX for some function A. This means that the potentials are 
stationary up to a gauge transformation. 

From = — i^dF + di^F = di^^F we also infer that there is a function (p so that 
i^F — —d(j). The function is a generalized electrostatic potential. Similarly, there is a 
generalized magnetostatic potential so that i^F — ~d(j), 

i^F^-d<j), i^F^-d(j). (18) 

It is obvious that these two generalized potentials are constant along the Killing trajectories, 
^^4> ~ and — 0. It is also clear that and are unique up to a constant. 

From the generalized Lorentz force equation — —CpiyF+igpiyF we obtain (Dyp) (^) = 
Dy{p{^)) — p{Dy^) — Cpiyi^F " igpiyi^F — —epv{(j)) + igpv{(l)). Since ^ is a Killing vector 
p{Dy^) = 0. This gives a conserved quantity 

— PiO + ^'p't' ~ i9p4' = const , (19) 

which is interpreted as total conserved energy consisting of a gravitationally modified ki- 
netic and rest energy p(^), a modified electrostatic energy Cpcf) and a modified magnetostatic 
energy igp4>. One may define a non-gravitational potential energy S'pot = ep4> — igp(f). Then 
— P{0 + <^pot- Note that the modified electrostatic and magnetostatic potentials are 
manifestly gauge invariant. 

The function cj) (the same holds for the function c^) is a generalized electrostatic (magne- 
tostatic) potential. Choosing an observer u collinear with the Killing field then dcf) = E 
where E = iuF is the electric field seen by the observer. From —dcf) = i^F = i^dA — 
£^(^A — di^A it is clear that in a gauge given by ^^A = we have (p = A{^). 

The Killing vector is defined up to a constant factor. The generalized electrostatic and 
magnetostatic potentials are defined only up to an additive constant. These two degrees of 
freedom allow the transformations 

^^^' = ae, 0->(/)' = (/)-)-/3, 0^0' = 0-Hz/3, (20) 

where a,/?, /? G IR with a 7^ 0. As a consequence, the total energy S" transforms as 

^ ^ ^' = -f a (ep/3 + gp$) . (21) 

Therefore, the factor a rescales the energy and a {cp/S + gpp) selects a new zero. This 
is in agreement with the experimental procedure of energy measurement where always 
differences of energies are compared - and such ratios of energy differences are in fact 
invariant under the affine transformations (|2ip . 

3.3 The model of the interferometer 

We assume that our interferometer is stationary and that the matter wave leaving the 
source is stationary, too. In more mathematical terms: 

1. The stationarity of the interferometer means that the worldline of each constituent of 
the interferometer is a Killing trajectory. Therefore, the 4-velocity of each constituent 
u is proportional to a Killing vector S^, u 

2. The stationarity of the source means that the kinetic momentum of the matter wave 
emitted by the source is time-independent with respect to the interferometer frame, 
^oPlsourcc ~ ^' where = denotes equality at the source in that particular frame. 

We now draw some consequences from these assumptions. 
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3.4 Consequences 

The condition of the stationarity of the source is equivalent to -S^cPlgourco ~ ^- Si^i^^s in 
a stationary space-time the covariant derivative and the Lie derivative commute, we have 
■^^P = everywhere. The latter also implies — 0. For the Lie-derivative of the 
canonical momentum we obtain = —CpdX + igpdX. This implies that the closed 

loop integral of the canonical momentum, the phase shift to be derived below, is time 
independent: 

^ j, ^ ^ ^ ^ -Cp ^ dX + igp ^ dX^O . (22) 

The momentum in the rest space of the interferometer is given by PuP- This 1-form on 
^ satisfies {PuP){£,) — and ^^P^p = 0. Therefore, P^p is also a 1-form in which 
we denote by p. The modulus of the momentum measured in ^ is |pp = h{p,p) [25]. It 
follows from the dispersion relation ([6]) 



\P\ = V {piOye-^U _ ^2 ^ ^(^ _ ^p^^)2 ^-2U _ ^2 _ (23) 

Since J^'^U = 0, -Sf^(/) = 0, and ^^4' = the functions U, (/), and (j) are also scalar functions 
on o/K. 

Note that uj^i, is a 2-form on the quotient space .yV and, thus, will be denoted by a;. 

4 Plebanski— Demianski space— time 

The initial form of the Plebahski-Demiahski metric in real coordinates describing the space- 
time with nonzero cosmological constant of a rotating and accelerating source endowed 
with mass, electric, magnetic and gravitomagnetic charges can be found in lllj . We follow 
the coordinate transformations presented in [571 dl] and write the metric of generalized 
accelerating black holes in the form 

ds^ = ^ (-^(dt-(asm'^9 + 2q(l-cos0))dip)^ + ^dr^ 
iV \ p'' A 

2 , 



with 



-L(^adt-ir^ + {a + qf)dipf +p^^^d9^^ , (24) 



^ + {q + acosey (25) 



a 



n = 1 (g + acos6')r (26) 



w 



P = sin^6'(l-a3COs6'-a4COS^6') (27) 

w V ui^ ' 3 



A = {K + e"^ + g^)-2mr + er^ -2n—r^ - [—K+-]r^ (28) 



where 



03 = 2a—m~Aaq^{K + e^+g^)-4:—aq (29) 

WW'' 3 

04 = -a'^^(n + e^ + g^)-^a^ (30) 

3 

e = -^^+Aq^m-{a^ + 3q')(^{n + e' + g^) + ^) (31) 
— g-^ w \w^ 6 ) 



— q'- w 
l + 2q^m-M%{e^+g^)~q^K 



(33) 
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We have 



.900 



A P 

— {asm^e + 2q{l~cos9)) - -^a{r^ + [a + qf) 



dip 



9ij 



-dr^ 



, sin^ 



1 



A 

(asin^9 + 2q(l~~cose)y 



^j,{r-^(^a^qff 



difi^ 



(34) 
(35) 

(36) 



This family of solutions are characterized by a mass-like parameter m, a cosmological 
constant A, a rotation-like parameter a, a NUT-like parameter q, the electric and magnetic 
charges e and g, an acceleration-like parameter a, and the twist parameter w. The latter 
can be given any convenient value provided both a and q are nonvanishing (all limits of 
this family of solutions including the limit of vanishing a or q and the corresponding choice 
of the parameter w can be found in |12|1. 

In order to obtain measurable components of physical quantities one has to introduce 
a tetrad defined by = u and CaCf^ = r/.j^, where T/.j^ = diag(— 1, 1, 1, 1). Henceforth tetrad 
components of a quantity are denoted by symbol ' . The components of the tetrad frame 
for the stationary observer are 



1 



* V^5oo 



/gr 
1 



1,0,0,0 
0,1,0,0 
0,0,1,0 



500/ 



/5;:;(^0, 1,0,0 
/5^fo, 0,1,0 



-900 



9av ~ 9009^^ 



gOy 

500 



,0,0,1 



/ gOy ~ gOOffyy 



'goo 



0,0,0,1 



(37) 
(38) 
(39) 

(40) 



Since all metrical coefficients do not depend on t, this metric admits a Killing vector 
= Therefore, 

-A + Pa2 



„2(7 



g(C,0 = goo 



and the acceleration of the Killing trajectories is 

1 



1 



a^, = d^jU = -df, ln(-5oo) = 7;9f, In 



'A-Pa2 



172p2 



(41) 



(42) 



Though the acceleration and the rotation of the Killing vector field can be calculated 
exactly, for the application to a realistic experimental situation on Earth the restriction to 
weak fields is sufficient and also leads to much simpler results. Therefore we assume the 
parameters m. A, a, q, a, w, e and g to be small and expand the metric coefficients to that 
order where these parameters appear first (e.g., first order in A and second order in e since 
there is no first order term for e). We obtain 

A . 1 



goo 



-1 



2m 
r 



{2q{q + a cos 9) - - g"^) - 2ra— cos 9 



and 



as 



m A , N 1 /„ , 9 9\ ct 
^ - — (r + to) + — (2q{q + acos9) - g^) + a— cos( 
3 r-' w 

aq sin C 



a — smt 

w 



(43) 



(44) 

(45) 
(46) 
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(It can be seen from (H^ that there are two constant terms in af- Both terms have recently 
been discussed as to whether they can be a possible origin of the Pioneer anomaly 28, 2S].) 

The rotation of the stationary congruence in the Plebahski-Demiahski space-time 
is given by 

^re = (47) 

'A m q 

3 r w ^ 

q f K m\ , „ ^, qa a 

+ ^ (g + 2a cos 6*) + -^ 

b r-^ I r w 



ujff, = a ( - - + ) sin6' (48) 
uj^^ = --^ + { ^ + ^] iq + 2acose) + —-COS0 . (49) 



If we regard the acceleration as the "electric" and the rotation as the "magnetic" part of 
the gravitational field, then the above expressions give rise to the interpretation that m is 
the gravitoelectric mass and q the gravitomagnetic mass. The product am may be regarded 
as a gravitomagnetic moment. 

Introducing the electric Afj_ and magnetic A^^ vector potentials for the metric (|24l) with 
components 

= (asin2 6' + 2^(1 - cos^)) - (g + a cos 6*) (r^ + {q + af) , 
p'^ ' ap^ 

. er + (q + a cos 6)g gr — (q + acos6)e 

At = 5 , At = I , 

p'' p^ 

K = (a sin^ 61 + 2^(1 - cos 61)) + ^{q^ acose){r^ + (q + af) , (50) 

we can write the total energy of the particle in the weak field approximation in the following 
form 

^ fc\ , fc\ , A - A fc\ , (^p^ + 9p9 , {cpg - gpe){q + a cose) 

S = p(0 + ^vot ^ P{0 + (ipAt - igpAt = p{0 + ^ + — ^^-^ , (51) 

where Cp and gp are electric and magnetic charges of the particle. 

For later use we note the measured components of the electromagnetic field 

Bf = ^ -2{q + acose) -- + 



2(q + a cos 9) (- 

(q + acosO) (-- - 4) , 





a 




a 




[r 


- q— 






w 




a 




a 






- q— 








w 



5 The phase shift 

5.1 The general expression 

The wave function propagates along two orbits, I and II, both starting at xo and ending 
at X, see also Fig. [2j The wave functions at those points are 

ipiix) = exp ( (feds) exp ( (f> ^ ) cq (54) 



2 

ipu{x) = exp ^-i j^6'(^s^ exp ^-i j^^^ flo, (55) 

where the upper integral is along path I and the lower along path II. The attenuation factor 
9ds will not play any role. The reason is that this factor only describes the attenuation 
of the density of the worldlines given by the phase S. However, no particle will disappear 
since we still have current conservation. At the analyzer the total intensity will be measured 
which is an integral over the probability density. This adds up to unity since nothing can 
disappear. Therefore we can omit the consideration of the attenuating factor. 
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The intensity of the two interfering wave functions at one port of the analyzer is 



1=^1+ = 2 (1 + cos A$) |ao|' , (56) 
where according to ([51)1 and ([55]) the phase shift A<i> is given by 

A$ = (j) ^ = j p + j [cpA - igpA) . (57) 

The second term is an Aharonov-Bohm hke phase shift and the third term its magnetic 
analogue. The first term gives the phase shift which stems from the change of the kinetic 
momentum due to the particle's interaction with the gravitational and electromagnetic 
field. 

We decompose the kinetic term p = p{u)u + P^p ~ p{^)e~^u + PuP and introduce the 
conserved energy (IT9l) 



= S" (j) e^^u ~ <j> [cpff) - igp(j)) e~^u + <j> PuP + <j> [cpA - igpA) . (58) 
Stoke's law gives 

A$ = S j e^^oj- j {^{cpdcj) - igpdip) A e^^ u + i^CpCj) ~ igpcj)) 

PuP + I {epF - igpF) . (59) 



With UHl) this yields 

A$ = ^ y e-'^Lo - j {ep(t) - igp$) e-^uj + j P^p + j P„ {epF ~ igpF) , (60) 

where we used the decomposition F = PuF + i^F A u (here P„ acts on both indices of F). 
Since PuP is a 1-form and lj as well as PuF and PuF are 2-forms on the quotient space 
jV we can write 

A$ = y e-^u; + j ^fpote-^o; + j> P + j {epF - igpF) . (61) 

where everything now is expressed in terms of the coordinates of . The last term describes 
the Aharonov-Bohm und dual Aharonov-Bohm effects [3D1I3I]. We introduce the magnetic 
and electric fluxes $m = J F and = i J F through the interferometer area. We also use 
and the unit 1-form n = p/\p\ in the direction of the momentum. Then the phase 
shift takes its final form 



A$ = / e"^a; + / (f. 



+ / \J (S" ~ S'potf e'^u - n + ep$™ - <?p$e , (62) 



opote 

It should be noted that this phase shift is manifestly gauge invariant and also invariant 
under the transformations (|20|) . 

The first two terms give a generalized Sagnac effect, the third term yields the phase 
shift due to the change of the potential energies over the interferometer, and the last two 
terms are Aharonov-Bohm like contributions. Here the Sagnac effect is modified by the 
gravitational potential and the non-gravitational potential energy. 



5.2 Small interferometer 

We now consider a small interferometer. Small interferometer means that all external fields 
like the gravitational potential as well as the electrostatic and magnetostatic potentials 
vary very slowly over the interferometer (Aharonov-Bohm like situations need a separate 
treatment). That means that the paths of the particles between the beam splitters, mirrors 



10 



mirror 1 




analyzer 



mirror 2 



beam 
splitter 



Figure 2: The gradient of U in the interferometer's rest frame, h and I are the interferometer's height and 
length. 



and analyzers are straight lines and the momenta of the particles are nearly constant. We 
also adjust the interferometer in such a way that the direction VJ7 is orthogonal to a line 
from the beam splitter to a mirror, see Fig. [21 Then path I is from the beam splitter — > 
mirror 1 — > analyzer and path II from the beam splitter mirror 2 — > analyzer. As far as 
the third integral in ((62)) is concerned we obtain 





/•mirroi^ \ 


^ mirrorl 


-/beam splitter / 



Ih ■ V J e-^u (^^ _ S'^ 



pot 



— '3pot,0 



beam splitter 



(63) 



where h is the vector connecting the beam splitter and the 1st mirror, / is the length of 
the interferometer and VC/ is evaluated at the position of the beam splitter (see Fig. 
The momentum of the particle at the beam splitter is po = [S — (?pot,o)^e~2'^<' — rrfl . We 
can choose the gravitational potential and the potential energy at the position of the beam 
splitter to vanish, [/q = and (S'pot,o = 0. Then 



-2(7 



i^pot)^ 



to2 n = -l—h-V(SV 

Po 



'pot; 



(64) 



Here and in the following we take the length I connecting the beam splitter and the 2nd 
mirror to be horizontal in the given coordinate system, that is, the beam splitter and mirror 
2 have the same r-coordinate. When a particle is not charged (fpot = 0. 

Non-relativistically, S" ^ m and po « mvo, where vq is the group velocity of the matter 
wave at the beam splitter. Then 



y^e-2^(^-^pot)2 



I 



-h-V{mU 



3pot 



(65) 



For vanishing Spot this is the result of [32] which has been confirmed by Colella, Overhauser, 
and Werner [5] using neutron interferometry. The appearance of the particle mass does 
not imply that the Weak Equivalence Principle is violated in interferometry or in quantum 
mechanics. Indeed, it has been shown by using quantum notions and measured expressions 
only that this phase shift can be given a natural form where the mass disappears. See |33j 
for a thorough discussion of this aspect. 

The first integral in ([62[) is the generalized Sagnac effect. It should be noted that it 
depends on the energy of the particle (frequency of the corresponding wave) only through 
the factor S". The integral itself is a purely geometric expression related to properties of 
the Killing congruence. For small interferometers we obtain 



(66) 
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where = ^e'^'^P'^u^ujpa- and = ^e^^p^u'^Y.P" (e'^"'"^ is the totally antisymmetric Levi- 
Civita symbol) are the rotation vector and area vector in .yV, respectively. 

Since we consider the case that the potential varies very slowly over the interferometer 
the potential energy can be put in front of the integral 



pote ^u; 



pot.O 



(67) 



so that for the case ^fpotjO = this term vanishes. (If we do not choose this condition, this 
term adds up with the first Sagnac term.) 

With these results we obtain as total phase shift 



A$ = ^ ( u; • S - — • V (<f [/ 

Po 



'pot. 



(68) 



where E and B are the electric and magnetic field strength given by Bi — Fjk and 
Ei = iFjk {i,j,k cyclic). This is the general result for a small interferometer. Now we 
have to insert the particular gravitational field given by the Plebanski-Demiahski metric 
and we also have to characterize the orientation of the interferometer. Since the baseline 
of the interferometer is always horizontal, there are only two degrees of freedom for the 
orientation of S. These two angles are the angle of the baseline with respect to the line 
9 = const, and the tilt of the interferometer. These two degrees of freedom then have to 
be matched to the longitudinal position of the laboratory on the surface of the Earth. 

Let 7 be the angle of the baseline with respect to and /3 the tilt angle. Then, in the 
local normalized 3-bein position of the laboratory 



h = h cos f3ef. — h cos 7 sin /Je^ — h sin 7 sin /3e^ 
S = E sin /3ef + S cos 7 cos /3eg + E sin 7 cos /3e^ , 

With this specification the phase shift ([68]) reads 



A$ = J'E 



(69) 
(70) 



— (cos Paf — cos 7 sin Pas — sin 7 sin Pa^) 
Po 



1 

Po 



[cos (idfS'pot - cos7sin/?9gf?pot - sin7sin/39^(?pot) 

+ sin Pujg^ + cos 7 cos Puj,^f + sin 7 cos /3w.g 

+epT,{smf3Bf + cos7Cos/3Bg + sin7C0s/3i?(^) 
— .gpE(sin jSEf + cos 7 cos f3Eg + sin 7 cos jSE^) , 



(71) 



where dfi — e'-^di,. This is the general result valid for a large class of interferometer orien- 
tations. 

From the fact that we can independently vary the angles /? and 7, we can extract from 
phase shift measurements the following combinations of terms 



A$(/3 = 


= 0,7 = 


0) 


= (fE 




S ' 
— af 
Po 


9r(opot ~ 

Po 


- CpTiBg — 


9p^Eg , 


(72) 


A$(/3 = 


f.7 = 


0) 


= S'Y. 




S 

PO \ 


H da^riot - 

Po 


- EpYiBf — 


gpYiEf , 


(73) 


A$(/3 = 


0,7 = 


f) 


= S'T. 




S ' 

Of 

Po 


dfS'pot 

Po 


- CpYiB^ - 


gpTiE^ , 


(74) 


A$(/? 


f>7 = 


f) 


= (fE 


^e^ + 


S 

— a A 

Po 


ffE 

H c/0(Spot 

Po 


■f CpTiBf - 


- gpTjEf ■ 


(75) 



Here the terms in the square brackets are gravito-inertial effects while the other terms 
are related to the electric and magnetic charges of the particles. In both cases we have 
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terms of similar structure: The first terms, e.g., af and dfS'pot are related to the change 
of the kinetic energy of the particle due to the interaction with the gravitational and 
electromagnetic field, while the other terms are Aharonov-Bohm type terms describing the 
coupling of the particle with the rotation flux or the magnetic and electric fluxes of the 
electromagnetic field. Note that each term is manifestly gauge invariant. 

It is obvious that charged particles have direct access to the electromagnetic field: an 
electrically charged particle has access to the magnetic field via the Aharonov-Bohm type 
terms and to the electric field via the derivative of the potential energy where in the latter 
the Killing vector has to be used. Both contributions can be separately identified through 
a variation of the momentum of the particles. The same holds, mutatis mutandis, for 
magnetically charged particles. For charged as well as neutral particles the contributions 
from the gravito-inertial field are of the same structure. In order to separate gravito-inertial 
and electromagnetic contributions one has to use two types of particles with different 
charge-to-mass ratios. 

Now we can use Eq. ([7T|) to explicitly calculate the phase shift for various observers in 
a Plebahski-Demiahski space-time. We start with an observer which is naturally aligned 
to the coordinate system used for expressing the Plebahski-Demiahski metric ([24|) . As a 
second example we take an observer who rotates in the given coordinate system. This can 
be accomplished by using a different Killing congruence. 



5.3 Adapted Killing observer 

The natural Killing observer is characterized by the Killing field = Sq in the coordinate 
system of (p4)) . From the previous subsection it is clear that with varying charges and 
velocities we have access to all components of aj, wjj, 9|^pot, E~^, and i?j separately. These 
quantities have been calculated in (|44| -(|49 | .([52 | . ((53|l . By further varying the latitude 9 
and the height r we then have independent access to the various parts of these quantities: 

m. A, Am, qa, 2q^ — e^~g^, a — , aq, aA, am, q, Aq, 

w 

O. (X Oi Oi 

mq, qa—, g, qg—, qe, ag—, ae, eaq—, e, qe—, qg , (76) 
www WW 

a a 
ae—, ag, gaq— . 

w w 

It should be noted that owing to the fact that the phase shift is a scalar calculated from a 
loop integral, it is independent of the chosen coordinate system. By varying the orientation 
and position of the interferometer and employing the expansion of the metric it is thus 
possible to extract all these parameter combinations. This particular procedure does not 
imply any interpretation of these parameters. In various physical settings other parameters 
and parameter combinations will appear 0. The fact that parameter combinations appear 
has been discussed to some extent in [121 El] in terms of deficit angles characterizing the 
conical singularities. Other parameter combinations appear if one calculates geometrically 
invariant parameters which, in general, do not have to coincide with the combination of 
parameters leading to physical effects like light deflection, perihelion shift, Lense-Thirring 
effect, etc. Our particular setup and the possibility to vary the experimental setting leads 
to the set of measurable parameter combinations (j76p . 

Combinations of these expressions give m. A, a, q, e, g, and a/w separately. If the 
experiment is carried out with neutral particles, we have access to and ljj- only and, 
thus, to the quantities in the first line of (|76p from which we can explicitly determine the 
parameter combinations 

TO, A, q, a, 5^ + , — . (77) 



^The same happens in, e.g., the discussion of effects within the PPN formaUsm: While for a 
Schwarzschild metric all effects are related to the mass (and this mass can also be given a purely geo- 
metrical meaning in terms of the surface of the horizon or a certain length-to— radius ratio) within the 
PPN formalism each physical effect like light bending, perihelion shift, red shift, etc, is characterized by a 
different combination of PPN parameters. 
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For neutral particles there is no separate access to e or 5. Also the acceleration a appears 
only through the combination with the twist parameter. 

The accuracy of the measurements of the acceleration is Aaf/flf w 10^^ [4J which 
on Earth implies an absolute accuracy of Aa « 10"^ m/s , and for rotation there is an 
absolute accuracy of buj « 10~^ s~^ [5]. In the near future the height on Earth and even 
to satellites orbiting the Earth at, say, 10000 km height can be defined with a precision 
of within 1 cm using Satellite Laser Ranging SLR f3F. This gives the relative error of 
distance estimations Ar/r « 10~^. 

All experimental results so far are in accordance with Newtonian gravity and Galilean 
kinematics both encoded in the mass parameter m and the rotation ui. That implies that 
the mass and the rotation of the Earth can be determined with an accuracy of 10~* and 
10""*, respectively. The accuracy is also good enough to detect gravitationally the presence 
of a mass of 100 kg located at a distance of ~ 1 m to the interferometer. 

Since all measurements are compatible with the Schwarzschild part of the Plebahski- 
Demianski metric we can place estimates, that is, maximum values, for the various Plebahski- 
Demianski parameters. Performing interference experiments on Earth as well as on a satel- 
lite 10000 km above the Earth surface and comparing the results would give an estimate 
A < 10^^^ m~^ (this is four orders of magnitude better than what one obtains from redshift 
experiments [MI)- If we take the NUT parameter q to be of dimension kg (which might be 
justified since it has the interpretation of a gravitomagnetic mass), then we obtain the es- 
timate q < 10^^ kg. From an asymptotic expansion of the metric and the requirement that 
the goi term should be of order we infer that a has the dimension m^/s. This leads to an 
estimate a < 10^^ m^/s~^. According to the previous settings the acceleration parameter 
a/w has the dimension m^^s^^ and, thus, we obtain the estimate a/w < 10""'^^ m^^s^^. 

5.4 Actively rotating observer 

Until now the interferometer is assumed to be positioned at a fixed point given by the 
coordinates (r, 9, (p) . In realistic situations the laboratory and, thus, the interferometer is 
attached to the surface of the Earth. Therefore both are co-rotating with the Earth. That 
introduces an additional angular velocity. This can be described as follows: Beside the 
time-like Killing vector ^ there is also a Killing vector related to the axial symmetry. In 
adapted coordinates this Killing vector is given by t]'^ =5^- The sum of two Killing vectors 
again is a Killing vector. Therefore, a stationary rotating laboratory on the surface of the 
Earth can be described by the Kilhng vector ^' = ^-l-w^Ty, where is the angular velocity 

^The gravitational field on Earth contains of course smaller additional contributions from the Sun and 
the Moon, for example. As a consequence one should consider the combined gravitational field leading to 
small modifications of the gravitational field considered here. However, since all the Plebafiski-Demiahski 
parameters except the mass are very small, the corresponding modifications of the gravitational field related 
to the parameters A, q, a, e, and g are safely negligible within our approximative scheme employed at this 
stage. Only modifications of the "Newtonian" part related to m have to be taken into account. These 
modifications can be calculated easily and are given, e.g., by ephemerides. It is no problem to include these 
modifications in our discussion of experiments. However, since we are more interested in the principle 
question of accessibility to the Plebariski-Demiahski parameters and in order not to blow up the formulas, 
we discuss the simpler question of measuring the gravitational field of the Earth only, neglecting all other 
masses in the Solar system. 
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of the rotating Earth. This new KiUing vector now gives a modified total rotation 

(78) 

- Ar (m+-) +2—r{2q + acose)] (79) 

\ 6/ W I 

l + -rA(2m + r)+r—(3q + acos0)] (80) 
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, (^"^ 








« sin 6* 1 
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(g + 2acos6l)-g-acos6' . (82) 
w ' 



as well as a modified acceleration 



2a;masin^ ( 



m rA 
^ " T 

- f r + m + -r'^AlSm + r) + 3—r'^a] sin^ 6* (83) 
\ 3 w J 

I ( 2m rA\ g 

1^^1/5 + w©asm26' — r- H — '2-^® sinfl 

« o y 3 / r 

/ A Q \ 

— ijjg sin Q cos ( r + 2m H — ^ (r + 4m) + 2-ar2 (84) 

a'^ = . (85) 

One easily recognizes the standard purely rotating observer in uj'^ , and io'f^ if one sets all 
parameters in the Plebahski-Demiahski solution to zero. 

This has to be inserted into the general result ((7T|) . As far as the access to the Plebahski- 
Demiahski parameters is concerned, we obtain the same result so that we omit to present 
the lengthy formulas. 



6 Conclusion 

We calculated the phase shift for a charged particle interference experiment in a general 
Plebahski-Demiahski black hole space-time. In doing so we put emphasis on a gauge 
invariant implementation of the symmetry conditions which are needed in order to obtain 
a stationary interference pattern. The gauge invariance ensures that each term in the phase 
shift is gauge invariant and, thus, has a clear physical interpretation. 

Besides addressing the issue of an appropriate formalism to describe such experiments 
we also answered the question whether it is possible to have access to all parameters charac- 
terizing the family of Plebahski-Demiahski generalized black hole solutions. The result ([77]) 
shows that the phase shift for a Mach-Zehnder interferometer in such Plebahski-Demiahski 
space-times is influenced by all parameters characterizing a general accelerating black hole 
defined by the metric ([24|) . By varying the orientation, the latitude and the height in the 
given gravitational field in principle it is possible to have access to all parameters of the 
Plebahski-Demiahski solution. 

The above results may also be used to describe atom interferometric experiments. In 
the non-relativistic limit the phase shift in a Kasevich-Chu like interferometer [36l |4] is 
given by A<I> = — fc • (a + t; x a;) where fe is the wave vector of the lasers serving as beam 
splitter and T is the time-of-flight of the atoms between two laser pulses, a and a; are the 
acceleration and rotation of the interferometer which can be taken to be the corresponding 
quantities calculated above. 
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